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Abstract. We present an algorithm to compute the GIT-fan of algebraic 
torus actions on affine varieties. 

1. Introduction 

Given an action of a connected reductive linear algebraic group H on an algebraic 
variety X, Mumford constructed open if -sets U C X admitting a good quotient 
U —> U//H, see [5D]. His construction depends on the choice of a iJ-linearized 
ample line bundle on X and, in general, one obtains several distinct quotients. 
This variation of GIT-quotients is described by a combinatorial structure, the so- 
called GIT-fan; see the work by Dolgachev/Hu [9] and Thaddeus 22! for ample 
bundles on a projective variety X and [3] for the affine case. 

In the present note, we provide an algorithm for computing the GIT-fan describ- 
ing the quotients arising from the possible linearizations of the trivial bundle for 
the case that X is affine and H is a torus. Note that the torus case is essential for 
the general one: if a connected reductive group G acts on X, then the associated 
GIT-fan equals that of the action of the torus G/G s on the affine variety X//G s , 
where G s C G is a maximal connected semisimple subgroup, see [5]. Moreover, our 
setting also occurs in the context of Mori dream spaces: there the Neron-Severi 
torus acts on the total coordinate space and the GIT-fan of this action is precisely 
the Mori chamber decomposition of the effective cone, see |14) . 

Our algorithm is based on the construction of the GIT-fan provided in [6j. We 
assume that X C W is given by concrete equations. The main computational 
steps are to determine the toric orbits of W intersecting X, see section [3j a suitable 
number of the so-called orbit cones of the i7-action on X and the GIT-chamber of a 
given weight, see section [2] The GIT-fan is then obtained by traversing a spanning 
tree of its dual graph; this idea also shows up in the computation of Grobner-fans, 
symmetric fans and tropical varieties as presented in [ITJ [TH [7] . We discuss some 
examples in section[4j At the moment, a Maple/convex [lOj implementation of our 
algorithm is available [T5] . 

The author would like to thank Jiirgen Hausen for valuable discussions and 
comments and the referee for helpful suggestions. 

2. Computing the GIT-fan 

Throughout the whole note, K is an algebraically closed field of characteristic 
zero. In this section, we first recall the necessary concepts from [6] and thereby fix 
our notation. Then we present and prove our algorithms for the GIT-fan. Aspects 
of efficiency of the algorithms are discussed at the end of this section. 

We will work with the following description of the toric orbits of K r in terms of 
faces of the orthant 7 := Q> : the standard torus T r := (K*) r acts via 

T x K r ->• K r , t ■ X = (tlXl, . . . ,t r x r ) . 
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Given a face 70 ^ 7, define the reduction of an r-tuple z, of e.g. numbers, along 70 
as 

._ ( / t \ t ._ e i € 70 

Z 7o - — 1*1) • ■ ' ) *W ) i " — 1 rv j 

[0, e.^70, 

where ei, . . . , e r € Q r denote the canonical basis vectors. Then, one has a bijection 

{ faces of 7 } <-> { TT-orbits } , 7o i-> T r Q := {i 7o ; i e T r } . 

Note, that in the notation of [12J . T 7o is the T r -orbit through the distinguished 
point corresponding to the dual face 7q := 7 1 n 7 V d 7 V - 

Definition 2.1. Let a C K[Ti, . . . , T r ] be an ideal. A face 70 of the positive orthant 
7 is an a-face if V(T 7o ; a) ^ 0. 

If X C K r is the zero set of the ideal a C K[T l7 . . . , T n ], then the a-faces cor- 
respond exactly to the T r -orbits intersecting X nontrivially. The computation of 
o- faces will be discussed in section OH 

We are ready to introduce GIT-chambers and the GIT-fan. Assume that the 
defining ideal a C K[Ti, . . . , T r ] oi X CW is monomial-free and homogeneous with 
respect to a Z fe -grading 

q t := dcg(T 4 ) eZ l , 1 < i < r . 

Then the corresponding action of the torus H = T k on W leaves the zero set 
X = V(W; a) C K r invariant. Let Q be the k x r matrix with columns q\,...,q r . 
We assume that the cone Q{j) C Q k is of dimension k. 

A projected a-face is a cone Q(jo) with 70 d 7 an a-face. In [6] these are called 
orbit cones. Write il a for the set of all projected o-faces. 

Definition 2.2. The GIT-chamber of a vector w e (5(7) = cone(gi, . . . ,q r ) C Q fe 
is the convex, polyhedral cone 

A(to):= p| 1JCQ 1 . 

u)Gi9Gfi„ 

The GIT-fan of the ii-action on X = V(K r ; a) is the set A(a,Q) = {A(to); to € 
Q(7)} of all GIT-chambers. 

As the name suggests, A(o, Q) is indeed a fan in Q k with (3(7) as its support, 
see [21 Thm. III. 1.2. 8]. Note, however, that the cones of the GIT fan need not be 
pointed in general. The set of j-dimensional cones of A(o, Q) will be denoted by 
A(a,Q)«. 

We turn to the computation of GIT-chambers. Let 51 := {Q(7o); 7o d 7} be the 
set of projected faces of 7 and let C fl be the subset of j-dimensional cones. 
Similarly, fi„ C il a is the subset of j-dimensional projected a-faces. We have 

n L a k) C fi^ := jtf € n (fe) ; all facets of 1? are in oi fe_1) | C , 

where the first containment is due to the fact that faces of projected a-faces are 
again projected a-faces, see [6j Cor. 2.4]. Given a vector to in the relative interior 
Q(7)°, set QW(to) for the collection of all 1? € that contain w. The next 

algorithm determines the associated GIT-chamber A = X(w). 
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Remark 2.3. (1) The set fi^') is computed directly by taking cones over suit- 
able subsets of {qi, . . . , q r }. 

(2) The computation of Q^(w) can be sped up via point location [18] . i.e. we 
only consider cones $ £ with at least one generator lying on the same 
side as w of a random hyperplane subdividing Q(j). 

(3) For an efficient computation of £2a , one reduces the amount of a-face tests 
as follows. Check for any d G £2™ if some 70 d 7 with Q(jo) = ?? is an 
a-face. As soon as such a face has been found, all other faces projecting to 
d may be ignored in subsequent tests. 

Algorithm 2.4 (GIT-chamber). Let w G Q(j)° be given. Assume that Q^ k \w) 
and £2 a are known. 

1 A := Q fc 

2 for each & G Q^(w) 

3 if 7? ^ A and all facets of d are in l ' 
5 return A 

Lemma 2.5. Let £ C Q fe 6e a pure k- dimensional fan with convex support |S| and 
let t £ I] be such that rfl |S|° 7^ 0. TTien r is i/ie intersection over all a G 5> fe ) 
satisfying t <a . 

Lemma 2.6. Let A G A(a,Q)W and z9 G £2[, fc) . n A° ^ fften A C tf . 

Proof. Suppose A % # - Choose it; G A° \ i? and u G n Then cone(u, w) n 
(i?o\^o) lies on some facet 770 ^ #o- By construction, ?7o n ^ 7^ 0- Since 770 G flif ^ 
holds, A is not a GIT-chamber; a contradiction. □ 

Proof of Algorithm \2.4\ The algorithm terminates with a cone A C Q k containing 
the given w G Q(j)° and our task is to show that A = X(w) holds. For this we 
establish 

A = 0^= fl ^ = A H- 



The first equality is due to the algorithm. The third one follows from Lemma 2.5 
Moreover, in the middle one, the inclusion "C" follows from Qfo' D Cl^ k '. Thus wt 
are left with verifying "D" of the middle equality. 

First suppose that X(w) is of full dimension. Then, for any i?q € ^0 with 



w G $0, we obtain i?g D X(w)° 7^ 0, because w G X(w)° holds. Lemma 2.6 shows 
X(w) C i?q. Thus, we obtain A 3 A(w). The case of dim(A(w)) < k then follows 
from the observation that X(w) is the intersection over all fulldimensional chambers 



X(w') with w G X(w'), see Lemma 2.5 □ 



Working with (fc — l)-dimensional projected a-faces in Algorithm 2.4 simplifies 
the necessary a-face tests compared to the following naive variant of the algorithm 
using fc-dimensional ones. 

Algorithm 2.7 (GIT-chamber, v2). Let w G Q("f)° be given and assume that 
^^(w) is known. 

1 A := Q k 

2 for each $ G ^^(w) 

3 if $ ^ A and there is an a-face 70 d 7 with Q(7o) = 1? 

4 A:=Afli5 

5 return A 
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The naive variant 2.7 in contrast, involves fewer convex geometric operations 
as |2.4| and thus can be more efficient if the latter ones are limiting the computation. 
See Remark |2.fl| for a more concrete comparison of complexity aspects. 

We turn to the GIT-fan. Given a full-dimensional cone A C Q fc , we denote by 
innerfacets(A) the set of all facets of A that intersect the relative interior Q{^f)° ■ 
Moreover, for two sets A, B, we shortly write A Q B for {A U B) \ {A H B). The 
following algorithm computes the set of maximal cones of the GIT-fan A(o, Q). 

Algorithm 2.8 (GIT-fan). 

1 A := {Ao} with a random full-dimensional GIT-chamber Ao 

2 J- :— innerfacets(Ao) 

3 while there is r\ G T 

4 Compute the full-dimensional GIT-chamber A' A with r\ -< A' 

5 A:=AU{A'} 

6 J- := J- Q innerfacets(A') 

7 return A 

Remark 2.9. In the fourth line of the algorithm, let A € A be the already found 



GIT-chamber with facet ?y. Then A' = A(u>') can be calculated with Algorithm 2.4 
where w' :— w(r]) — e ■ v for some w(i]) E rj° and u € A v fl i) 1 with a suitably small 
e > 0. One possibly must reduce e until X(w') n A = r\. 

Proof of Algorithm \2.8[ Write |A| for the union over all A € A and \T\ for the union 
over all r\ <E T . Then, in each passage of the loop, a full-dimensional chamber of 
A(o, Q) is added to A and, after adapting, |J r |nQ(7)° is the boundary of |A|nQ(7)° 
with respect to Q(j)°. The set T is empty if and only if |A| equals Q{~f). This 
shows that the algorithm terminates with the collection of maximal cones of A(a, Q) 
as output. □ 



Note that Algorithm 2.8 traverses a spanning tree of the (implicitly known) 
dual graph of A(a, Q) which has the maximal cones as its vertices and any two are 
connected by an edge if they share a common facet. Another traversal method for 
implicitly known graphs is reverse search by Avis and Fukuda [3] , which also might 
be applied to our problem by the following observation. 

Proposition 2.10. The GIT-fan A(a,Q) is the normal fan of a polyhedron. If 
Q>o ^= Qil)) then k{a,Q)( k } can be enumerated using reverse search. 

Proof. The first statement is [TJ Cor. 10.4]. The second claim follows from the first 
one and [TTJ Sec. 3]. □ 



Remark 2.11. We compare the usage of Algorithm 2.4 (in 2.8 ) to that of 2.7 As a 



test, we compute the GIT-fans of the maximal torus action on the (affinc cones over 
the) Grassmannians G(2,5) and G(2, 6), using a Maple/convex implementation. 
The following table lists the total number of a-face tests and the total number of 



cones $ entering the fourth line of Algorithms 2.4 and 2.7 



G(2,5) 
G(2,6) 



Algorithm 
J a-face-tests 
30F 
6574 



2.8 



with 2.4 

() cones •& 
21 
50 



Algorithm 2. 
(j a-face-tcsts 
469 
21012 



2.7 



with 

j) cones 
2(T 
52 



Note that in Algorithm |2.4| the a-face tests concern faces of lower dimension than 
in Algorithm |2.7| and thus are even faster. 



Remark 2.12. (1) Intermediate storage of occurring cones and their intersec- 

saves time. 



tions in Algorithms 2.4 and 2.7 



(2) The traversal of the GIT-fan can take advantage of symmetries as explained 
in PI Ch. 3.1]. 
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3. Computing o-faces 

Let X C W be the zero set of an ideal a C K[Ti, . . . , T r ]. Here we compute torus 
orbits of K r intersecting X. In the notation of section [2] this means to determine 
the a-faces 70 d 7 of the orthant 7 = Q> . 

Given a face 70 d 7 and a polynomial / g K[Ti , . . . , T r ] , we write / 7o := f(T~ ) 6 
K[T 7o ] where T :— (Ti, . . . , T r ), i.e. we replace each Tj with zero if ^ 70. Let 
a 7o '■— (/70) / *= a ) — ^[^70]- A direct a-face test is the following, based on a 
radical membership problem. 

Remark 3.1. A face 70 < 7 is an a-face if and only if Y[ e £7o T± g" ^<x l0 . 

This leads to a Grobner based way to decide whether a given 70 -< 7 is an a- 
face. The main aim of this section is speed up this direct approach by dividing 
out all possible torus symmetry. This is done in A lgorithm |3.5| Further possible 
improvements are discussed at the end of the section. 

First consider any torus T and an ideal c C 0(T). Let H C T be the maximal 
subgroup leaving V(T; c) invariant and denote by 7r: T — > T/H the quotient map. 
To describe tt explicitly, we use the correspondence between integral matrices and 
homomorphisms of algebraic tori: every n x k matrix A defines a homomorphism 
a : T k -> T" by sending t g T k to (t Al * , . . . , t A " ) g T" where the A it are the rows 
of A. 

Remark 3.2. The map tt: T k = T — > T/H = T" is given by any n x k matrix P 
of full rank satisfying 

ker(P) - f) ker ( P s)' 

gee 

where to g = a^T 110 + . . . + a m T Vm € c we assign the m x k matrix P g with rows 

Remark 3.3. Fix a generating set G := (31, ...,ji)ofcC K[T±\ T* 1 }. Let 
Pg be the stack matrix, i.e. the vertical concatenation, of P gi , . . . , P g! . Compute 
the Hermite normal form D = U ■ Pq with a unimodular matrix {/. Choose P as 
the matrix consisting of the upper non-zero rows of D. Then P describes tt: T k = 
T -> T/H ^ T". 

Proof. Clearly, P is of full rank. Since the exponent vectors of each g £ c are linear 
combinations of the exponent vectors of g% , . . . , gi, we have 

l 

ker(P) = ker(P G ) - f| ker(P ffi ) = f| ker(P s ) . 

i=l gGc 

□ 

A push forward of g € c under tt is a h € 0(T/H) satisfying n*h = T^g for some 
monomial T M ; we simply write ir^g for any such /i and 

tt*c:= (tt^; .9 g c) C G(T/P) . 

Remark 3.4. Let P: Z fc -> Z" be as in [O and let g = a a T"« + ... + a m T v ™ g c. 



B := P g ■ V ■ a ■ U 1 where er is a rational section for D and P g is as in 3.2 Then 



Compute a Smith normal form D = U-P-V with unimodular matrices U, V. Define 
B := P g ■ V ■ a ■ U- 1 where 
there is /1 G Z> such that 

tt^ = T M (oo + aiT Sl - + . . . + a m T s ™') g K[T lt . . . ,T n ] . 

Proof. Let tt 9 : T k -> T m and /3 : T n -> T m be the maps of tori defined by the 
matrices P g and B. Clearly, g = T K TT* g h for h := a + a x Ti + ... + a m T m g 0(T m ) 
and some k g Z> . Each 3 g c is P-homogeneous. This implies ker(P) C ker(P 9 ), 
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so there is a unique integral matrix B' such that P g = B' ■ P. In particular, B = B' 
is integral. Therefore, g = T K ir*(/3*h). □ 

We now specialize to the case of a-face-verification. Given 70 di 7, let #(70) C 
TT be the maximal subgroup leaving V(T^ o ; a 7o ) invariant. Our approach reduces 
the dimension of the problem by using 

V(T r lQ ; a) ± & V (TjH( 7o ); tt,o 70 ) ^ 0. 

Algorithm 3.5 (o-face verification). Let o = (fi, . . . , f s ) C K[2i, . . . , T r \ be an 
ideal and let 70 -< 7. Set gi := (/i) 7o and G :— (gi, . . . , <7 S ). 

1 Use|3. 3|to compute a matrix P representing 7r: T r — > T 7o /iJ(7o) 

2 Apply |3.4| to P to obtain n*G := (zr* 51 , . . . , 7r*<7 s ) 

3 ifTi---T„ e yplCK^,..,!,,] 

4 return false 

5 return TRUE 

Proof. The map n is a good quotient for the iJ(7o)-action on T 7q . Consequently, 
we have 

\i=l / i=l 

by standard properties of good quotients [HI p. 96]. This shows that V(T^ o ; a 7o ) 7^ 
if and only if V(T n ; tt^G)) ^ 0. □ 

Remark 3.6. If the total number of terms occurring among the generators is 
low compared to the number of variables in the sense that P = Pq in the first 
line of Algorithm |3.5[ then we might speed up the algorithm using linear algebra 
as follows. Each term n*gi is linear by construction. Solve the linear system of 
equations 71% G = 0. Then 70 is an a- face if and only if there is a solution in T™. 

Let us briefly recall the connection to tropical geometry, compare e.g. [7|. Given 
a monomial-free ideal a C K[Ti, . . . , T r ], its tropical variety is 

trop(a) := f]trop(f) C Q r , 

where trop(/) is the support of the codimension one skeleton of the normal fan of 
the Newton polytope of /. By |21) . 

(1) 70 d 7 is an a-face trop(a) H (70) ^ . 

Fixing a fan structure on trop(o), this can be turned into a computable crite- 
rion. Note however that trop(a) usually carries more information than needed to 
determine the a- faces and is in general harder to compute (see [7] for an algorithm). 

Remark 3.7. To compute all a-faces, the number of calls to Algorithm |3.5| can be 
reduced by any of the following ideas. 

(1) The tropical prevariety of a generating set (/1, . . . , f s ) of o is the coarsest 
common refinement n^Ti where is the one-codimensional skeleton of the 
normal fan of the Newton polytope of fi. Then each face 70 ;_ 7 whose 
dual face 7g does not satisfy equation w.r.t rijYj is not an a-face. 

(2) A face 70 ^ 7 is not an a-face if and only if there is / € a such that exactly 
one vertex of the newton polytope of / lies in 70; also compare [SJ Prop. 
9.3]. Choosing any subset of a, we may identify some faces 70 d 7 that are 
no a-faces. 



COMPUTING THE GIT-FAN 



7 



(3) Veronese embedding: Let 70 ri 7 be such that there are (classically) homo- 
geneous generators gi, . . .,g 8 of o 7o of degree d e Z> . The images of the 
gi under 

K[T 70 ] -> K[S M ; Ml + • • • + Mr = <*] , T M H> 5 M 

give a linear system of equations with coefficient matrix A. If a Gauss- 
Jordan normal form of A contains a row with exactly one non-zero entry, 
70 is no o-face. Adding redundant generators to a 7o refines this procedure. 

(4) Let <t e S r be a permutation of (the indices of) the variables Ti,...,T r 
that keeps the set of generators of a invariant. Then 

7o d 7 a-face cone(e CT (i); € 70) o-face . 

Some of those permutations can be computed by assigning a both edge- 
and vertex-colored graph to the generators of a and computing its auto- 
morphism group, e.g. using |19j . 

Remark 3.8. The efficiency of Algorithm |3.5| depends on the algorithms used 
for both Grobner bases and Smith normal forms. An implementation using the 
respective built in functions of Maple gave the following timings. 



remark 


3.1 


algorithm 


3.5 


with 


3.7 


;2) 


a-faces of d2,5 
a-faces of 02,6 
a-faces of 02,7 
a-faces of 02,3,3 


21s 
16 min 
> 3 days 
4.03 h 


10s 
76 s 
24.8 h 
44.1 min 







There, 02, n stands for the respective Pliicker ideal and 02,3,3 denotes the defining 
ideal of the Cox ring of the space X(2, 3, 3) of complete rank two collineations [T3l 
Thm. 1]. 



4. Examples 



We consider torus actions on the affine cone over the Grassmannian G(2, n) 
induced by a diagonal action on the Pliicker coordinate space K r , where r — Q). 
Such actions will be encoded by assigning the variable Tj the i-th column of a 
matrix Q. Moreover, we write 02, „ C K[Ti , . . . , T r ] for the Pliicker ideal. 

We compute both, the GIT-fan of the torus action on V(K r ; 02, „) as well as the 
GIT fan of the ambient space W . The latter coincides with the so-called Gelfand 
Kapranov Zelevinsky decomposition GKZ(Q), i.e. the coarsest common refinement 
of all normal fans having their rays among the cones over the columns of Q. In 
general, the Gelfand Kapranov Zelevinsky decomposition is a refinement of the 
GIT-fan. See [5] for a toric background. 

Below, the drawings show (projections of) the intersections of the respective fans 
with the standard simplex. 



Example 4.1. For n = 4, the ideal 02,4 
homogeneous with respect to 



(T.Te - T 2 T 5 + T 3 r 4 ) CK[T U ...,T ( 



6 IS 



Q 





A(o 2 ,4,Q) GKZ(Q) 

Using Algorithm |2.8[ we obtain the four maximal GIT-chambers of A(o2,4,Q). 
The finer fan GKZ(Q) has twelve maximal cones. 



Example 4.2. 

respect to 



For n = 5, the ideal 02,5 C K[T 1; . . . , T 10 ] is homogeneous with 
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Q = 



By Algorithm 



1 
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1 
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1 
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1 
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1 





1 











1 








1 





1 


1 





A(a 2 , 6 ,<9) GKZ(Q) 

there are twelve four-dimensional cones in A(a2,5, Q) whereas 



GKZ(Q) contains 336 such cones. 

Example 4.3. For n = 6, the ideal d2,6 ^ K[Ti, . . . , T15] is homogeneous with 
respect to 



Q 



1 














1 
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1 




















1 











1 











1 
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1 


1 
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1 














1 











1 








1 





1 1 



Using Algorithm 2.8 we obtain the 81 five-dimensional cones of K{a2fi,Q)- The 



fan GKZ(Q) has 61920 such cones. 
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